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Abstract 

In this paper we give a new Koksma-Hlawka type inequality for Quasi- 
Monte Carlo (QMC) integration. QMC integration of a function f: [0,1)® —>■ 
R by a finite point set V C [0,1)® is the approximation of the integral 
^(/) ■= /[o.i)'= by the average IvU) ~ Exsp /W- We treat 

a certain class of point sets V called digital nets. A Koksma-Hlawka type 
inequality is an inequality bounding the integration error Err(/;P) : = 
Hf) - Ivif) by a bound of the form |Err(/;P)| < C ■ ||/|| • D{r). We 
can obtain a Koksma-Hlawka type inequality by estimating bounds on 
|/(k)|, where /(k) is a generalized Fourier coefficient with respect to the 
Walsh system. In this paper we prove bounds on Walsh coefficients /(k) 
by introducing an operator called ‘dyadic difference’ di^n- By converting 
dyadic differences di^n to derivatives g|-, we get a new bound on |/(k)| for 
a function / whose mixed partial derivatives up to order a in each vari¬ 
able are continuous. This new bound is smaller than the known bound 
on |/(k)| under some condition. The new Koksma-Hlawka inequality is 
derived using this new bound on the Walsh coefficients. 


1 Introduction and the main results 

Quasi-Monte Carlo(QMC) integration of a function /: [0,1)® ^ K by a finite 
point set V C [0,1)® is the approximation of the integral /(/) := /(x) dx 

by the average Ivif) '■= 1^ Ex6-p/(^) O’ CH] and [22] for details). We 
want to find quadrature point sets V making the absolute value of the integration 
error |Err(/;T’)| := [/(/) — Iv{f)\ small for a set of functions /. This problem 
is formulated as follows: We consider a function space H with a norm ||/||_fr 
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and the worst case error su]^||^<i |Err(/;7’)| by a QMC rule using the point 
set V (for example, see [9], [l^ for details). Then, it holds that, for any f € H, 

|Err(/;lP)|<||/||ff X sup |Err(/;iP)|. (1) 

ll/lln<l 

Thus in order to make |Err(/;P)| small, we have to obtain quadrature point 
sets V making the worst case error sup|jj||^<;^ |Err(/;T’)| small. 

We often treat a point set V called ‘digital net’ (for example, see [H]). A 
digital net V is defined as follows. Let n, m, 6 > 1 be integers with n > m. Let 
0 < /i < 6™ be an integer and Ci,... ,Cs he nxm matrices over the finite group 
Zf, = Z/6Z. We write the 6-adic expansion h = take a vector 

h = [hi ,..., hm) G (Z™)^, where hj is considered to be an element in Z;,. For 
1 < i < s, we define the vector {yh,i,i, • ■ •, yh,i,n) = h • (Ci)^ and a real number 
Xi{h) = J2i<j<n yh,i,jb~^ € [0,1), where yh,i,j is considered to be an element of 
{0, ...,&— f} ^ Z. Then we define a digital net V by {xq, • • • , where 

x/t = We define the dual net [SI H], which is essential to 

analyze the integration error: 

7?^ := {k = (fci,..., ks) GWo\Cjki + --- + Cjh = OG Z™}, 

where ki = (k^p, ..., for ki with &-adic expression ki = J2j>i ■ 

Here Kij is considered to be an element of Zf,. Throughout this paper, when 
we take a point set V, we assume that 7^ is a digital net with 6 = 2. 

In the classical theory, many researchers studied the integration error of a 
function / with bounded variation (or function with square integrable partial 
derivatives up to first order in each variable) (for example, see |9], [16]). An 
extension to smooth periodic functions was established in [3], while a further 
extension to smooth (non-periodic) functions was shown in [4]. The QMC rules 
constructed in these papers, called higher order QMC rules, achieve (up to pow¬ 
ers of log N) the optimal rate of convergence. See also [5] for more background 
on higher order QMC rules. The purpose of this paper is to substantially im¬ 
prove the constants in the bounds on integration error in [3], which is crucial 
in problems in uncertainty quantification [Zldl]. In particular, m point out 
that the large constants from [4] cause problems in the CBC construction of 
interlaced polynomial lattice rules, which is one of construction methods to ob¬ 
tain point sets whose worst case error achieves the optimal order. To avoid this 
problem, they suggest to use much smaller constants which are more realistic. 
This paper provides the theoretical justification for doing so. 

We explain the details. Dick et al. [7] introduced a smooth function space 
whose functions / satisfy that their norms (|4|) (see below) are finite. If / is 
a function whose mixed partial derivatives up to order a in each variable are 
continuous, then / is contained in this space. This space has some parameters 
called weights C M>o = {x G K : x > 0} ,where S' := {1,..., s}, which 

model the importance of different coordinate projections, see [23] . 

To state their results, we need modified dual spaces which correspond to the 
subsets V C S. For k.;, G let (kt,; 0) denote the vector whose jth component 
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is kj if j € V and 0 otherwise. We define the dual space which corresponds to 
the subset (j> d S hy Vy := {k„ G | k = (k„; 0) G V^] (note that none 
of the components in is 0). Let 1 < r, r', g < oo with l/r + 1/r' = 1 and 


Ma((^lj • ■ • J ^|t)|)) — ^ ^ 1 ) 

2 = 1 j<OL 

for li with dyadic expansion li = with > ■•• > ai^Nf They 

showed the following bound on the worst case error (Dick et al. [7] also showed 
the results for a digital net with b > 2): 


sup |Err(/;T’)| < es^a,-fy('P), 


r<l 


with 


/ 


Ijr' 


Cs,a,7,r''(^) — 


k;„G-PG 


This implies the following inequality of the form (HD: 

|Err(/;T’)| < \\f\\s,a,j,q,r X es,a,j,r'{'P), 

where 

||/IU.™.:= (e E 

\uCS \ iiCu T„\„G{l,...,a —1}I“\"I 


/ / dys\v 


r/g\ 


( 2 ) 

(3) 

(4) 


1 /r 


with the obvious modifications if g or r is infinite. 0 Here (a„, Ty\y , 0) denotes 
a sequence {vj)j with Vj = a for j G v, Vj = tj for j G u\u, and Vj = 0 for 
j ^ u. And we write /("i. ■."<») = a’"i+-+’"<*//aa:7 • • • 

Based on these bounds on the integration error, Dick constructed ‘interlaced 
digital nets’ to obtain a point set with small integration error (for example, see 
MM)- He showed that the worst case error of this type of point set achieves 
the order 0(Af““(log Ny°‘) in terms of the cardinality A^ of a point set (see [4]). 
This is known to be optimal up to log terms (see |21)1. In mm there is also a 
component-by-component (CBC) algorithm to obtain point sets which achieve 
the same order. 

There is another algorithm to find good quadrature point sets for QMC 
for integrands with large enough smoothness a. This was introduced by Mat- 
sumoto, Saito and Matoba m- They define the Walsh Figure of Merit (WAFOM), 

^The norm in [3 Definition 3.3] has been corrected in arXiv: 1309.4624 ^3. The correct 
version is restated here in Eq. Q. 
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which is defined by the discretization of the upper bound of the worst case error 
of the form ©■ The advantage of WAFOM is that we can compute it on the 
computer in reasonable time. This property enables us to find a point set with 
small integration error by computer search. In fact, there are some algorithms 
for finding good point sets for QMC (see [13],[14]). 

In this way, to find good quadrature point sets, we need an inequality of 
the form o, which bounds the integration error by the product of a norm of 
/ and a figure of merit of V. These types of inequalities are called Koksma- 
Hlawka inequalities (for example, see [16] for details). In the following, we give 
a new Koksma-Hlawka type inequality to bound the integration error of smooth 
functions better than the inequality ([3]) under some condition. 

Theorem 1.1. Let a G NU{oo} such that a >2. We assume that a function f 
satisfies that its mixed partial derivatives up to order a in each variable Xi are 
continuous on [0,1]®, and 1 < P-,q,q' < oo such that l/g + 1/q' = 1. Then we 
have 


|Err(/;iP)| < x W„.^,,(iP), 


where 


q\ 1/9 


Wa.,,,(7’) =1 E h- E 

y<j>^vcs \ k„e'P,t 


1/9' 


(^«^2V sup 


and where 




dMf 

— - dyLs\y 

.ui ax, 




P \ i/p 
dx„ 


with the obvious modifications if either p, q or q' is infinite. 
In this theorem we write 

|i;| min(Q;,A^i) 

Ma((^ij • ■ ■ I ^|«|)) = E E (“hi + 2) 

i=i i=i 


(5) 


for li = 2“''^ instead of Dick’s weight function HaHh, ■ ■ ■, /|u|))- 

This result yields a significant improvement of ([3]). This is crucial when 
using the bound in a CBC algorithm, since a large constant (as it appears in 
[3 Theorem 3.5]) may make it impractical to perform the CBC construction in 
practice. For instance, m Section 4.1] write that The resulting large values of 





Bounds on Walsh coefRcients 


5 


the worst-case error bounds [ referring to the large constants in [71 Theorem 3.5] 
] have been found to lead to generating vectors with bad projections. 

Additionally, we also include the case a = oo which has not been studied 
before in the context of digital nets. In [71 Theorem 3.5], the case a = oo 
is not included since in this case the constant Ca appearing in ([3]) is infinite. 
Furthermore, we can define another version of WAFOM when we consider this 
new bound (see m for details). 

This theorem is based on the estimation of Err(/;7^) by the Walsh coeffi¬ 
cients. Dyadic Walsh coefficients are defined as follows (see [ID],[20] for details). 

Definition 1.2 (Walsh functions and Walsh coefficients). Let f: [0,1)'* — )• K. 
and k = (ki,... ,ks) G Ng. We define the h-th dyadic Walsh function walk by 

S 

walk(x) := J^(—l)(^j>i 

i=l 

where for 1 < i < s, we write the dyadic expansion of ki by ki = 

and Xi by Xi = > where for each i, infinitely many digits bij are 0. 

Using Walsh functions, we define the h-th dyadic Walsh coefficient /(k) as 
follows: 


/(k) := / /(x) • walk(x)dx. 

d[0,l)'* 

We see that the integration error Err(/; V) by a digital net V can be repre¬ 
sented by Walsh coefficients /(k) as follows ([H Chapter 15]): 

Err(/;7^)= ^ /(k) (= E E />C 0 )). 

keP-L\{o} 

The proof of Theorem 11.11 is facilitated by an improved bound on the Walsh 
coefficients of smooth functions. We show bounds on Walsh coefficients /(k) as 
follows: 

Theorem 1.3. We assume the same assumptions as in Theorem \l.l[ Let (j) yf 
u C {1,..., s} = A. For k^ e we have 


|/((k„;0))l < 2^ • ||/(min(a,Nk„))| 


P) 


( 6 ) 


where 1 < p < oo and jj • jjp is the norm defined in Theorem \l.l[ Here we define 
the symbol min(a, Ni) = (min(a, A^i),..., min(a, Af|«|)) for 1 = (/i,..., Z|„|) with 
dyadic expansion f = 2“*’" ■ 

This inequality follows from the formula for the Walsh coefficients by dyadic 
differences, which are defined in Section [3| (see the rough sketch of the proof in 
Section [3]) . 
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Here we compare this result with [SI Theorem 14] and its higher dimensional 
analogue in |7] . Note that our bound includes the case a = oo for the case 5 = 2 
and we see that our bound m is better under some condition. Assume that 
s = 1. Then for A^i > a, if we multiply our bound by (5/3)““^ our bound is still 
smaller than that bound by Dick for any ki = (see [U chapter 14]). 

If A^i < a it is in general not clear which bound is better. 

In the following, we assume that A: S Nq has dyadic expansion with k = 
2“^ where N is some integer and ai > • • • > aw and set = 0, {uj} = (j) 
for A: = 0. 

The remainder of this paper is organized as follows. In Section |2 we give 
the proof of Theorem [IT] and we give the rough sketch of the proof of Theorem 
11.31 in Section [Sj In Section |4| we show the proof of lemmas to complete the 
proof of Theorem 11.11 and 11.31 


2 Proof of Theorem 11.11 

Proof. We assume that / is continuous on [0,1]® and 
[H Lemma 17], we have pointwise absolute convergence 


/(x) = ./(k)walk(x). (7) 

keN* 


Now we have that / is continuous by the assumption of /. We show that / also 
satisfies the second condition. If we apply Theorem 11.31 for a = 2, we have 


E !/>)!= E E l/((k.;0))|< E E 2^2-^^(*^”)||/»n.J) 

kGNg\{0} 0/t;CSk„eNl’’l k„eNl”l 


< 


(p^vCS k„GNl’'l 


max 

nG{l,2}l*’l 


11/^^ 


|]„ < 2p max 

n'g{0,l,2}® 


11/^’^ 


\LP 


E 

4>^v<zS k. 




GNI”! 


Note that L^-norm 11 /|]lp := (/[o i)<i 1/(2^) is different from the norm 

II/Up defined in Section 1. Thus we have 


E i/>)i 

kGNg 


= l/(0)l+ E !/>)! 

kGNg\{0} 

< 2 ” 


max ||/("')||ip(l+ V V 
n'G{0,l,2}'» V ^ ^ / 


b/«CSk„GNl’'l 
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Since niaXn/g{o,i, 2 }'’ ll/*’" ^IIlp < oo holds by the assumption on /, we have only 
to show the last summation is finite. We prove this in the following way: 

1-f ^ ^ = A + 

k^eNl”! V fceN / 


= h + E2 


-1-2 4- ^ ^ 


= 


-Z 2-4 


\ 

Then we can apply the formula © to / to get 

1 


|Err(/;lP)| = 


/(x) dx - ^ /(> 

bo,!)* I' I 


xeP 


< 77 + 2 


1 


-4 


nf! 






/(O) - T^T E E /(k)walk(x) 

' ' xer kGN;; 


Now we introduce the property of Walsh coefficients walk- Let "P be a digital 
net in [0,1)® where |P| = 2™. Then we have (see [9j Lemma 4.75]) 


walk(x) = 

xev 

Using this fact, we have 


2^ if k e r^, 
0 otherwise. 


|Err(/;P)| = 


/(o) - E />) 

xGP-L 


< y] i/(k)i= y] y] i/((k„;o))|. 

xG^-^Vto} <l>¥^vcS 


Let a S N U { 00 } with a > 2 and 1 < p,q,q' < 00 such that 1/q + l/q' = 1. 
Applying Theorem 11.31 to /, we have 


|Err(/;P)| < y] y] 2^2-^:^(*^'') ||/(“' 


(min(a,Nk„))| 


^ k„G'P,f 
M 


<E7„^2p sup ||/^“•’^||p 

y a„e{l,...,a}l"l k„e'P-L 


— II/IIBq, 7.P.<?' ^ LVa,7,ij(P). 

We use Holder’s inequality in the last inequality. 


□ 


3 Proof of Theorem 2 

3.1 Some notations 

To calculate bounds on Walsh coefficients, we introduce the dyadic difference 
di^n and the weight function /ru(k). 


< 00 . 
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Definition 3.1 (dyadic difference). Let s,n,i G N with i < s. For a function 
g: [0,1)® R, we define the dyadic difference di^n{g) by 


di,n{g){xi, ...,Xs) 


g{xi, ...,Xi®2 - g{xi, ..., Xj,..., x^) 

2-n 


Here we write z © 2“" := z + 2“"(—1)^" for z having dyadic expansion z = 
> where infinitely many digits Zj are 0. 

Let a vector k = {ki,... ,ks) G Ng, u = (ui,..., Ug) G (No U {oo})® with 
ki = 2“'’^ ■ Let (f) V G S = and a vector k„ G with 

k = (k^,; 0). The symbol (k„; 0) is the same symbol as in Section 1. We use the 
symbol dk„ instead of the composition map of {^i^aiJ+l}i^v,l<j<min{Ni,Ui}■ 

'Eiemar]^ 3.2. Since any two dyadic differences commute, dk„ is defined inde¬ 
pendent of the order of a composition. 

Definition 3.3 (the new weight function /i^(k)). We use the same symbols as 
in Definition \3.1[ The weight function /ru(k) o/k is defined by 


Mu(^) ^ ^ ff 2), 

i^v,l<j<min.{Ni ,Ui} 


and we define /ru(0) = 0 . 

When iti = a G N U {oo} for every i, ^^((kt,; 0)) equals p'^fky). 

We define the important two functions Xn{x,y) and W(k). 

Definition 3.4. We define the function Xnix, y ): [0,1)^ —>• R 6j/ 

, , f 2"’ if y G [min(a::, x © 2“"), max(a::, x © 2“"’)1, 

Xn(x,y):=| Q otherwise. 

Recall that x © 2“" is defined as in DeRnition \S.l[ 

Using this, we define the 1-dimensional versions W{k) inductively by 

W{0) := 1, 

W(2"i)(y) := [ Xn^+iix,y) dx, 

Jo 

W{2^- + • • • + 2"-+i)(y) := [ Xn.^,+i{x, y)W(2"^ + • ■ • + 2^-){x) dx, 

Jq 

where ni > ■ ■ ■ > n^+i. 

Then, for a vector (fci,..., kg) G Ng, the s-dimensional versions W(k): [0,1)® -G 
R is defined by 


W(k) -Jlwih). 
2=1 
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Remark 3.5. By definition, W{k) is continuous on [0,1) for any k € Nq. 

The following important property of these functions is proven in Section [4.3l 

Lemma 3.6. Let a vector fc G N. We have that W{k) >0 on [0,1) and, for 
1 < p < oo, we have 

\\W{k)\\Lv <2^-^. 


3.2 Proof of Theorem 11.31 

In this subsection, we show bounds on the Walsh coefficients /(k) admitting 
the following Lemmas 13.81 and 13.91 which we prove in the next section. 

We denote the symbols used in the statements. 


Definition 3.7. We use the same symbols as in DeRnition \3.1[ For k and u, 
we define 


i,> 


T.j>u, 2“^'^ if i&v, 

0 if i & S\v, 


2“^'^ if i&v, 

0 if i & S\v, 


and 


k" := := 

ki," := k^,< := 

min(u,Nk) := (min(ui, min(u,Nk„) := (min(iti, 

I min(u,Nk)|zi := y^/min{ui, Nj). 

iGS 


When we analyze Walsh coefficients, it is suitable to use dyadic differences. 
In fact, the k-th Walsh coefficient /(k) can be represented by 9i^„/(k) as follows: 

/(k) = (-1) • 2-“---2a.,a7^(/)(k), (8) 

where k = {ki... ,ks) with ki — 2“*’^. Applying the above formula ([5]) 

repeatedly, we have the following formula: 

Lemma 3.8. Let f G L^([0,1)'*), k G Ng and u G (No U {oo})'*. 

Then we have dk„/ G L^([0,1)®) and we have 

/(k) = (_l)l min(u,Nk)|,l (9) 

This formula (|9]) means that dyadic differences connect the k-th Walsh co¬ 
efficient /(k) to the weight function ^^(k) for / G I)'*)- 

Dyadic differences 5^ „/ are similar to derivatives In the following for¬ 
mula, we can replace di^nf with using a slight modification. 
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Lemma 3.9. Let u = (ui,... ,Us) G (Nq U {cx)})'*. We assume that a function 
f satisfies that its mixed partial derivatives up to order Ui in each variable Xi 
are continuous on [0,1]®. For any vector k G Nq, we have 

/(k) = (_l)|min(u,Nk)|,l 2 -M:(k) f /('"“("■Nk)) (j,) . ) (x) ■ walfc- (x) dx. 

Then we can get the following bound on |/(k)|: 

Lemma 3.10. We assume the same assumptions as in Lemma 1,9. ,91 Let (j) ^ 
uC{l,...,s} = 5' and k^ G NI’'! . Then we have 

|/((k„;0))| < 2^ • < oo, 

where 1 < p < oo and || • ||p is the norm defined in Theorem \l.l[ 

Proof. We write x„ = {xi)i^v for x G [0,1)®. We see that W((ki,; 0)" )(x) = 
W(ki,<)(x„) and wal(k„;o)“ (x) = walk„u(x„) since W(0) = walo = 1. )And we 
see |walk„“ | = 1 by the definition of Walsh functions. Combining those facts 
and Lemma 13.91 we have 


|/((k„;0))| 


r ^((min(u,N,J;0))(^) . ) (x) ■ walk„^ (x) dx 

Jlopy 


'[ 0 . 1 )” 


'[0,1)IS\H 


/((m”(u,NkJ;0))(^) dxs\. 


|w(k^<)(x„)| dx.„ 


< 2-ML((k.;0))||j(min(u.NkJ)||^ . ||VL(k^li; 


Lp-1 


where we used Holder’s inequality in the third inequality. By Lemma 13.61 we 
have 


||W(k 4 )|| 


Thus we obtain the result. 


□ 


In particular, when rti = a G N U {oo} for every i, we have Theorem 11.31 
In the following section, we will prove the lemmas which we used in this 
section. From now, we denote by 0"=! Pi composition of maps tpi o • • • o 


4 Proof of Lemmas 

Definition 4.1. We use the same symbols as in Definition \3.1\ . Let p = 
(Pi)iei>,q = {qi)iev G with I < Pi < Qi < Ni We use the following symbols 
in the proof. 

n di^aij+1- 

iG.v,pi<j<qi 

We use this symbol when we can recognize k. 
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4.1 Proof of Lemma 13.81 

Proof. We prove only the case s = 1 here. In the case s > 1, we obtain the 
result by applying the same method in a component-wise fashion. 

We easily obtain the first statement as follows. Let ki = 2“^ • Since 

di^aj+i is the sum of f{xi ©2““^“^) e L^([0,1)) and / S L^([0,1)), we have 
di,aj+if G L^([0,1)'*). By repeating this argument, we have ^(/) e 

L\\Q^)). 

We show the second statement inductively. We omit the case fci = 0 or 
ui = 0 since the proof is easy. We show the case mi = 1: 

/(fci) = (-1) • (10) 

where ki = 2“^ ■ By changing variables cci cci © 2““^“^, we have 


+ 


[ /(xi©2 ^) ■ walk^ixi) dxi 

Jo 

G(2c+1) 

f{xi +2““^“^) ■ walfei(xi) dxi 
■2c 

/(xi -2““^“^) •wal/ci(xi) dxi^ 

G(2c+2) 

/(xi) • walfei (xi - dxi 

<2c+l) 

/(xi) • walfei (xi + 2““^“^) dxi'j 

= f /(xi) • walfci(xi © 2““^“^) dxi = f /(xi) ■ walfei(xi) ■ walfei(2““^“^ 
Jo Jo 

= - f{xi) •waUi(xi) dxi, 

Jo 

where the last two identities follow from the definition of Walsh functions. Using 
this calculation, we obtain 


-1 

E 

c^O 

(E 

r 

J"G(2c+2) 

J2““J 

-G(2c-H1) 

2“J -1 

E 

c^O 

(E; 

r 

J ^■(2c+l) 

J2~‘‘3 

“^■2c 


di,aj+i{f) (ki) = -2 ■ 2“^+^ • J /(xi) • walfci (xi) dxi = (-1) • 2“-^+^ • /(fci). 

We write U = min(ui, W)- Using (|T(I1) inductively, we obtain 

(fci) = d[f)V (*i) = (-1) • 2“^+' • dgV (fci) 

= (- 1)2 • {h) = ■■■ = (- 1 )^ • 


which is the result. 


□ 
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4.2 Proof of Lemma 13.91 

4.2.1 Important properties of dyadic differences di^n 

In order to prove Lemma 13.91 we show some properties of dyadic differences 
9i,n- We define the following symbols. 

Definition 4.2. We use the same symbols as in DeHnition \3.1\ and \4.1\ We 
define 

Wi^n := wal2n-i(a;i), 

wdi^nig) ■■= Wi^n ■ d^^n{g) for g [0,1)® -)> M, 

and 

:= wd^ := 

Notice that is a function but wdp is an operator. We can rewrite the Walsh 
function as follows: walk = 

We see that wym and di^n commute in the next lemma. 

Lemma 4.3. When {i, n) {j, m) G for a function g : [0,1)® —>■ M, we have 
the following identity: 


‘^j,m ‘ ^i,n{g') — di^n{g ' 

Proof. We omit the proof here since it is easy. □ 

We first prove the following property. 

Lemma 4.4. We use the same symbols in the above definition. For a function 
g G L^([0,1)®), we have wdpg G L^([0,1)®). 

Proof. Since the above lemma, we have that wd^t/ equals • d^^. By the 
definition of Walsh functions, we see |w^| = 1. And since d^g is the sum of the 
functions in L^([0,1)®) as in the proof of Lemma [3l8l we have d^g G L^([0,1)®). 
Thus the result follows. □ 

4.2.2 Proof of Lemma 13.91 

The following Lemma is the key to prove Lemma |3.91 which connects wdi^nig) 
with the derivative 

Lemma 4.5. Let n, s,i € N satisfy s > i. Let g: [0,1]® —>■ R as a function of 
the ith component Xi, satisfy 

n — 1 2 


g € C^[2-^+^c,2-^+\c+ 1))) 


c = 0,...,2’ 


( 11 ) 
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Then for any z = (zi,..., Zg) G [0,1)®, we have 
f ^ ^9 * 

= J {zi,...,z^-i,y,z,+i,...,Zs) ■Xnizi,y) dy, 


where we define 
dxi dxi 


on X 


i G [2-”+ic, 2-"+i(c +1)), c = 0,..., 2"-i - 1. 


Proof. Let c' G No satisfying Zi G [2“"c', 2“”(c' + 1)). We consider two cases: 
c' = 2c and c' = 2c + 1 for some integer c. We only calculate the case c = 2c' 
since the other case can be calculated by the same way. In this case, by the 
calculation Wi^n{zi) = 1 and the assumption CD, we have 


wdi^n{g){'z) 


= 9i,n{g){z) = 2” • {g{zi,.. 




■, Zi + 2 , ■ ■ ■, Zg) g{zi,..., Zg)) 

Zi-i,y,Zi+i, ...,Zg)dy 



Zi-l,y,z^+l,...,Zg) -Xnizr^v) dy. 


The last equality follows from [z,, Zi+2 "] = [min(zi, z^ © 2 ”), max(zi, Zj © 2 ")] 
and the definition of □ 

We prove Lemma [T9] using these results. 


Proof. We prove the case s = 1. We omit the case fci = 0 or ui = 0 since the 
proof is easy. We write fci = 2“-’ and U = min(ui,iV) here. We assume 

that ui > 1, then we have 

0 

We use Lemma 14.31 in the third equality. Using the assumption of / and the 
definition of Wa^+i, we have that 


/ • G ([2-“''+i-ic,2-“^+i-i(c+l))), 


and we have 


— ( 


( df 


dxi V‘ 

f ■ ^(^U+l)) 

H 

II 

• W(G+1)J 


with 0 < c < 2 “a+i+i _ I 

Let 1 < n, 0 < c' < 2” be integers. By the definition of wdi^n, we have that, 
if g G ([02“", (c+ 1)2“")), it holds that wdi^ng G ([c2“", (c+1)2“")) and 
-£^{wdi^ng) = wdi^ni-^) on [c2“", (c+ 1)2“"). 
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If we take g = f ■ yv 


(N) 

iU+1) 


and n = 0(7 + 1, we have 


— on [c2 S{c+1)2 i), 

where 0 < c < — 1. Applying this argument inductively, we have 

^ = wdj 27 (^ on [c2-“^-\(c+l)2-“=-i), 


where 0 < c < 2“^+^ — 1. Since 2 ^ > 2 we can take n = oi + 1 and 

9 = (/ ■ ■''^(( 7 + 1 )) in Lemma [4.51 Then we continue the computation of 

follows 


0 

Now we have S L^([0,1)) since |rya^+i| = 1 and the assumption of 

/. Therefore if we take g = in Lemma [4.41 and consider the fact 

|Xai+i(a;i,y)| < 2“i+\ we see that the integrand Xai+i • (^■£^ ■ 

in the last line belongs to L^([0,1)^). Thus we can use Fubini’s Theorem as 
follows. 

Xai+i(a;i,y)dxi^ • 

= • W(aii))(y) ■W{2‘^^){y) dy 

= ■w[|^|^))(a:i)-IT(2“^)(xi)dxi. 
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By repeating the argument we have 

(ki) = wdi,a 2 +i (a;i) • W^(2“i)(a;i)da:i 

= (/ ^“2+i(a;i>y)-W^(2“')(a:i)c^a:i^ 

= / ^d [37 (g . wW) {y) ■ W^(2- + 2-)(z/) dj; 

^ lo +2“")(a;i)da:i 

^ Jo {'JJ 

thus we have the result for s = 1. By calculating in a component-wise manner, 
we have the result for the case s > 1. We omit that case here. □ 

In fact, we can determine the sign of /(k) in the special case. 

Corollary 4.6. Let f G C°°([0,1]®) and k G Nq. We use the symbol Ni appear¬ 
ing in Definition \3.1[ Then, have /(k) • (—l)5Df=i'^i > Q. 

Proof. By Lemma ld^ and the fact W(0) = 1, we have W (k) = rifci ^(^0 0. 

By combining this fact and the assumption that -Wd > .^g have that 

the product . \y{\P) > 0. Thus, by Lemma 13.91 with Ui = oo we have 

/(k) • (-l)^?=i^* = • f /W.....JVd(x) . W(k)(x) dx > 0, 


which is the result. □ 

4.3 Proof of Lemma 13.61 

4.3.1 Important properties of Xn(x,y) and W{k) 

We show the important properties of Xn{x,y) and W{k) in this subsection. See 
Definition 13.41 for the definitions of Xn{x,y) and W{k). 

We see that Xn{x,y)/2'^ is a characteristic function of some region in [0,1)^. 

Lemma 4.7. Let n S N and c S Nq satisfying c < 2”“^. 


1. Let x,y G [0,2 then we have 


Xn{x + c2 ”+\y-tc2 "+i) = x„(a;,y). 
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2. Let X € [c2-"+i, (c + 1)2-^+^) and y ^ [c2-^+\ (c + l)2-’"+i). Then we 
have Xn{x,y) = 0. 

And let y £ [c2-^+^, {c + 1)2-^+^) and x ^ [c2-^+\ {c + 1)2-^+^). Then 
we have Xn{x,y) = 0 . 

Proof. 1. We have (a: + c2“"+^)©2“” = (a:©2“") + c2“"+^. Thus, the result 
follows from the fact that 

y £ [min(a:, x © 2“"), max(a:, x © 2“")] 

y + c2-”+i e [min(cc + c2-”+\ {x + c2-^+^) © 2“"), 
max(x + c2-’^+\ {x + c2-"+i) © 2-”)]. 


2. We prove Xn{x,y) = 0 in the case x £ [c2“"’+^, (c + 1)2“"+^) and y ^ 
[c2“"+^,(c + 1)2“"+^). Let X £ [d2“", (d + 1)2“") for d £ Nq where 
d = 2c or 2c + 1. When d = 2c, it holds that c2“"+^ < x < x © 2“"’ = 
X + 2“" < (c + 1)2“"’+^. In the case d = 2c + 1, it holds that c2“”+^ < 
X © 2“" = X — 2“" < X < (c + 1)2“"+^. So we have 

[min(x,x©2-”),max(x,x©2-")] C [c2-”+i, (c + l)2-”+i). 

So, if y ^ [c2“”+^, (c+1)2“”+^), we have y ^ [min(x, x©2“”), max(x, x© 
2“”)]. Then we obtain Xn(x,y) = 0. 

For the case x ^ [c2-”+i, (c + l)2-"+i) and y £ [c2-”+\ (c + 1)2-”+^), 
there is some integer e such that x £ [e2“"+^,(e + 1)2“"+^) and y ^ 
[e2“"+^, (e + 1)2“"+^). Thus the result follows from the above argument. 
We finish the proof. □ 

The function Xn(x,y) is defined by using a characteristic function of y. In 
Lemma [T8l we rewrite Xn(x,y) by using a characteristic function of x. 

Lemma 4.8. Let y £ [0, 1) and c, n S Nq satisfy y £ [2“"c, 2“”(c+ 1)). 

If c = 2c' for some integer c', we have 


Xn{x,y) 


2" if xe [2-’"c,y]U[2-"(c+l),y + 2-”], 
0 otherwise. 


And if c = 2c' + 1 for some integer c', we have 


Xn{x,y) 


2" if X G [y - 2-", 2-"c) U [y, 2-"(c + I)), 
0 otherwise. 


Proof. We only prove the case c = 2c' here since the case c = 2c' + 1 follows 
from the same argument. In this case, by Item 2 of Lemma 14.71 we have that 
for all y G p-^+^c', 2-"+i(c' + 1)), 

Xnix, y) = 0, X ^ [2""+ic', 2“"+i(c' + 1)). 


( 12 ) 
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Let X £ [0,1) and d £ No satisfy x G [2“"d, 2~'^{d + 1)). We calculate Xn{x, y) 
in the three cases: d ^ {2c', 2c' + 1}, d = 2c' and d = 2c' + 1. 

We consider the case d 0 {2c', 2c' + 1}. By condition (|T^ . we have 


Xn{x,y)=0, x£[2 "d,2 ”(d+l)). 


In the case d = 2c', since x®2 "=x + 2 we have 


Xn(.x,y) = 2" 


X <y < x + 2 ”, 

2 -n+ic' <x,y< 2-”(2c' + 1). 


So we have 


Xn{x,y) 


2” if X e [2-”+ic',y], 

0 if xG (j/,2-”(2c' + l)). 


When d = 2c' + 1, by a similar argument to the case d = 2c', we have 

, , r 2” if xG [2-”(2c' + l),?/ + 2-”], 

Xn{x,y) if a; g (j; + 2-”,2-”(2c' + 2)). 


By combining these cases, we have the result. □ 

In the last lemma, we show the period of a function W{k). 

Lemma 4.9. We have that W{k) is a periodic function with period 2~°‘^ for a 
positive integer k = 2“'- 

Proof. We proceed by induction on N. We prove the result for k = 2“L Let 
c G No satisfying c < 2“L We have that for y £ [0, 2““i), 


W(2“i)(y + c2 “1) = [ Xai+i{x,y + c2 “i)dx 
do 

= / Xa,+i{x,y + c2-°-^)dx= Xa^+l{z + c2-°■\y + c2-°■^) 

Jc2-^i Jo 

= f Xa^+iiz,y)dz= [ Xai+i(2,y)dz = W(2“^)(y). 
do do 


dz 


The second and fifth equalities follow from Item 2 of Lemma 14.71 the forth 
equality follows from Item 1 of Lemma 14.71 and the change of variable x = 
z + c2~‘^^ in the third equality. 

Now we assume that the lemma holds for the case k]\{ = Sti We prove 
the result for the case k = satisfying qn+i < un. By the induction 

assumption, we have that W{ki\f){z + d2~°'^) = W{ki\f){z) for z G [0,2““") and 
an integer d satisfying 0 < d < 2“". Let c G No satisfying c < 2“"+L Then we 
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have that for y g[0,2 


Wik){y + c2-“'^+i) = WikN + 2“^+^)(?/ + 02““^+^) 

= [ XaN+i+i{x^y + c2~°'^+^) ■ W{kN){x) dx 
Jo 

^(c+l)2““'V+l 

XaN+i+i{x,y + c2"“"+i) ■W{kN){x)dx 


/ c 2-^ n+i 
2-“Af + l 


Xa^+^+iiz + c2-“«+\ J/ + c2-“«+i) • Wik^Kz + c2-“^+^) dz. 


The third equality follows from Item 2 of Lemma 14.71 and the change of variables 
x = z + 02““"^+! in the last equality. By the induction assumption, we have 
W{k]sf){z + c2““"+i) = W{kN)iz) for z G [0, 2““"+i). Thus we obtain 


lT(fc)(y+ c2-“'^+i) 


2-“n-+i 


XaM+i + l{z + C2- 


■O-N + 1 


y + c2 ■ W{kN){z) dz. 


Then we continue the computation as follows: 


2-“iV+i 

W{k){y + c2-‘^^+^)= / Xai^+i+i{z,y) ■W{kN){z)dz 

Jo 

= [ Xa!,+,+i{z,y) ■W{kN){z)dz = W{k){y). 

Jo 


The first equality follows from Item I of Lemma 14.71 and the second equality 
follows from Item 2 of Lemma 14.71 □ 


4.3.2 Proof of Lemma 13.61 

We prove Lemma 13.61 using the results in the above subsection. 

Proof. Since Xn > 0, we see that W{k) >0 on [0,1) by induction. We omit the 
details. We use this property of W{k) >0 to prove ||II^(fc)||Lp < 2^^“p^. 

Using Holder’s inequality, we have 

\\Wik)\\%= [ \Wik){x)\-\W{k){xr-^\dx<\\W{k)\\L4Wik)\\l^^^ 

J[0,1) 

Thus we have 


wwmLP <\\wik)\\l\\wik)\\lJ. 

Then, if we have ||IU(fc)||ii = 1 and ||IT(/c)||l=o = 2, we have 

\\Wik)U. < \\W{k)\\l\\W{k)\\]-J =2^^-^^ 
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Therefore we prove ||lT(A:)||ii = 1 and ||lT(fc)||Loo = 2 for any A: G N to complete 
the proof. 

We prove the case k = by induction on N. In the case fc = 2“i, we 

have 


||W(2“^)|Ui 


n Xai+iix,y) dydx 


fo JO 


0 J 


2°-^+^ dydx = 1. 


The first equality follows from Xai+i{x,y) > 0. We prove ||W(2 “i)||l~ = 2 to 
complete the case k = 2“^. By Lemma 14.91 we have 


||W(2“0IIl~= sup \W{2-^){y)\= sup \W{2-^){y)\. 

yG[0,l) !/g[0.2-“i) 


Since W(2“i)(?/) = Xai+i(x, y) dx and Xai+i(x,y) > 0, we have 

||IL(2“i)||l- = sup f Xai+i{x,y)dx 
i;eI0.2-“i) Jo 


= max 


sup / Xar+i{x,y)dx, sup / Xai+i(a;,y) 

\!/G[0.2-“i-1) Jq i/G[2-“i-i,2-“i) Jo 


dx . 


We calculate the supremum on [0,2 and [2 ^,2 “i) separately. We 

assume that y G [0,2““^“^). By Lemma IT^ we have 

/•l ry ^2““i“^+y 

/ Xai+iix,y)dx = / 2“i+^da;+ / 2“^+^ dx. 

Jo Jo J2-“1-1 


>0 Jo J 2 - 


If we choose y = 2 we can maximize the right hand side. Hence we obtain 


sup 

ye[0,2-“i-i) 



Xai+i{x,y) dx 



= 2 . 


By the same argument we get the same result in the case y G [2 ^, 2 . 

We omit the details. Therefore we have ||IT(2“i)||ioo = 2. 

Now we assume that for any 2“S we have that ||IH(fcjv)||j;^i = 1 

and ||IT(A:Ar)||Loo = 2. Let k = k^ + 2“"+^ satisfying qn > qn+i- We prove 
that for any k = k^ + 2°'^+'^, l|bI^(^)llLi = 1 and ||IT(fc)||Loo = 2. By the fact 
that WiltN + 2“'^+'^) > 0 and Fubini’s Theorem, we have 


||W(/cw + 2“"+i)|Ui 



XaN+i+i{x,y)W{kN){x)dxdy = 

Jo 


W{kN){x) dx. 


By the fact W{kN) > 0 and the assumption on k^, we obtain ||IF(A:Ar + 
= \\WikN)\\L^ = 1. 
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We prove ||W(fc)||j;,oo = 2 as follows. By the fact W{kN + 2°-^+^) > 0 and 
Lemma [4.91 we have 

||W(fc)||Loo = sup W(fcw + 2“^+^) 
yg[0.2"“w+i) 

= max ( sup / XaN+i+i{x,y)W{kN){x)dx, 
^ye[0,2-“w+i-i) Jo 

sup 

yg[2-“w+i-L2““«‘+i) 

We can calculate the supremum by the fact W{k]sr) > 0 and the same method 
as in the case k = 2“^: 

suPyg[o, 2 -“^^+i-i) lo Xa„+i+i(a;, y)WikN){x) dx 
supj^g[ 2 -<x^,+i-i 2 -;v+i) fo Xa^+,+i(x, y)W(k;v)(x) dx 

Then we obtain 


2-«iV+i 


2^^+i+W(A:Ar)(x) dx. 


XaAT + i + l 


{x,y)W{kN){x)dx'^ . 


2-«iV+i 


\mk)\u 


2“"+i+W(fcjv)(x) dx = 2“"+ 


2“W-“N+1_1 2-“n 


1 + 1 


2“iV “N+1_1 


E 

2=0 


^ -j. iV « 

2 a;,+i+i / W{kN)ix)dx = 2^^+^ ■ 

i=o "'o -^0 


22““^ 

2-“iV 


TT(A:Ar)(^) c?x 
W{kN)ix) dx. 


The second equality follows from Lemma 14.91 Thus we have 


pZ J'' 

||W(fc)||ioo =2“"+i- / WikN){x)dx 
Jo 

W(kN){x) dx = 2 • f W(kN){x)dx. 

Jo 


2“^-l .(2+1)2““^ 

= 2 E [ 


The first equality follows from Lemma |4.9I By the assumption on few and the 
fact W{kpf) > 0, it follows that W(kj\i)(x) dx = ||W(fcAr)||Li = 1, and hence 
we obtain ||W(fc)||L=o = 2 • W(kjY)(x) dx = 2. □ 
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